Probability and complex quantum trajectories 



Money V. John 

Department of Physics, St. Thomas College, Kozhencherry, Kerala 68964-1, 

India. 



Abstract 

It is shown that in the complex trajectory representation of quantum mechanics, 
the Born's ^>*^> probability density can be obtained from the imaginary part of the 
velocity field of particles on the real axis. Extending this probability axiom to the 
complex plane, we first attempt to find a probability density by solving an appropri- 
ate conservation equation. The characteristic curves of this conservation equation 
are found to be the same as the complex paths of particles in the new representa- 
tion. The boundary condition in this case is that the extended probability density 
should agree with the quantum probability rule along the real line. For the simple, 
time-independent, one-dimensional problems worked out here, we find that a con- 
served probability density can be derived from the velocity field of particles, except 
in regions where the trajectories were previously suspected to be nonviable. An al- 
ternative method to find this probability density in terms of a trajectory integral, 
which is easier to implement on a computer and useful for single particle solutions, 
is also presented. Most importantly, we show, by using the complex extension of 
Schrodinger equation, that the desired conservation equation can be derived from 
this definition of probability density. 

Key words: quantum Hamilton-Jacobi equation, trajectory representation, 
probability axiom, complex methods 
PACS: 03.65.Ca 



1 Introduction 



In a previous work [T], the quantum Hamilton-Jacobi equation [2-5] (QHJE) 
was made use of to demonstrate the existence of particle trajectories in a 
complex space, for different quantum states. This complex quantum trajec- 
tory representation was obtained by modifying the de Broglie-Bohm (dBB) 
approach to quantum mechanics [6], which allows particle motion guided by 
the wave function. One of the advantages of the resulting theory, which offers a 
new interpretation of quantum mechanics, is that it does not face the problem 
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of stationarity of particles in bound states, encountered in the dBB represen- 
tation. Another trajectory approach to quantum mechanics, which also claims 
the absence of this problem, is the representation developed by Floyd, Faraggi, 
Matone (FFM) and others [7-9]. 

The new complex trajectory representation proceeds by first substituting \& = 
e iS/h j n Schrodinger equation to obtain the quantum Hamilton- Jacobi 
equation 
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and then postulating an equation of motion 
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for the particle. The trajectories x(t) of the particle in the complex x-plane are 
obtained by integrating this equation with respect to time p] . It was observed 
that the above identification \1/ = e tS l h helps to utilize all the information 
contained in \1/ while obtaining the trajectory. (The dBB approach, which uses 
^ = Re tS l h does not have this advantage.) The complex eigentrajectories in 
the free particle, harmonic oscillator and potential step problems and complex 
trajectories for a wave packet solution were obtained in [JJ. The representation 
was extended to three dimensional problems, such as the hydrogen atom, by 
Yang [10] and was used to investigate one dimensional scattering problems and 
bound state problems by Chou and Wyatt [HH2]- Later, a complex trajectory 
approach for solving the QHJE was developed by Tannor and co-workers |13j . 
The QHJE was derived independently by Sanz and Miret-Artes [J4] , who also 
found the complex trajectory representation useful in better understanding 
the nonlocality in quantum mechanics pT5TfT6] . 

It is well known that the QHJE as given in Eq. was used by many physicists 
such as Wentzel, Pauli and Dirac, even during the time of inception of quantum 
mechanics [2]. In a commendable work in 1982, Leacock and Padget [4] have 
used the QHJE to obtain eigenvalues in many bound state problems, without 
actually having to solve the corresponding Schrodinger equation. However, 
there were no trajectories in their work and it was only in [1] that the equa- 
tion of motion (TS]) explicitly solved and the complex trajectories of particles 
in any quantum state obtained and drawn, for the first time. This paper also 
highlighted the interpretational value of the complex quantum trajectory rep- 
resentation, vis-a-vis the Bohmian mechanics. Eventhough this formulation 
requires that the wavefunction is known from a separate calculation, it is gen- 
erally conceded that it was the work in [I] which provided a complex trajectory 



2 



interpretation of quantum mechanics p^6]fT7] . 



Instead of computing the complex trajectories x(t), the complex paths Xi(x r ) 
in the above scheme can directly be found by integrating the equation 

dXi X{ . 

d/Xi- Xr> 



where Eq. ([2]) shall be used. In [T| , it was noted that even for an eigenstate, the 
particle can be in any one of its infinitely many possible quantum trajectories, 
depending on its initial position in the complex plane. Therefore, the expec- 
tation values of dynamical variables are to be evaluated over an ensemble of 
particles in all possible trajectories. It was postulated that the average of a 
dynamical variable O can be obtained using the measure as 

oo 

<0>= [ OV*V dx, (4) 



where the integral is to be taken along the real axis PQ. Also it was noted 
that in this form, there is no need to make the conventional operator replace- 
ments. The above postulate is equivalent to the Born's probability axiom for 
observables such as position, momentum, energy, etc., and one can show that 
< O > coincides with the corresponding quantum mechanical expectation val- 
ues. This makes the new scheme equivalent to standard quantum mechanics 
when averages of dynamical variables are computed. 

One of the challenges before this complex quantum trajectory representation, 
which is an ontological theory of particle motion, is to explain the quantum 
probability axiom. In the dBB approach, there were several attempts to obtain 
the probability distribution from more fundamental assumptions [6]. In 
the present paper, we first attempt to obtain this distribution along the real 
line from the velocity of particles in the complex trajectory representation. 
It is found that always there exists a direct relationship between the 
distribution and the imaginary component of the particle's velocity on the 
real line. Since this distribution is defined and used only along the real axis, 
the conservation equation for probability in the standard quantum mechanics 
is valid here also, without any modifications. 

At the same time, since we have the complex paths, it would be natural to 
consider the probability for the particle to be in a particular path. In addition, 
we may consider the probability to find the particle around different points in 
the same path, which can also be different. Thus it is desirable to extend the 
probability axiom to the x r Xj-plane. But in this case, it becomes necessary to 
see whether probability conservation holds everywhere in the plane. A recent 
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paper by Poirier [18] addresses this issue and obtains some negative results for 
the choices made for such a distribution. Poirier tries to define an extended 
probability by working with the solution of time-dependent Schrodinger 
equation. The probability density at a point x in the complex plane is pos- 
tulated as p(x) = ty(x)^(x) where \fr(x) = \fr*(x*). The complexified flux is 
chosen as j(x,t) = v(x,t)p(x,t) = —(ih/m)ty*(x*)ty' where v(x,t) = x is 
given by equation (j2J) and prime denotes spatial differentiation. With the help 
of time-dependent Schrodinger equation, the author shows that, in general, 



This, arguably leads to nonconservation of probability along trajectories. But 
it shall be reminded that this negative result is based on the choices made in 
[T8] for the probability density and flux. 



In the present paper, on the other hand, we first show that the Born's prob- 
ability density, defined along obtainable as the exponential of an 
integral over the real line, of Xi, the imaginary part of the particle's complex 
velocity. Next, extending this probability axiom to the entire x^-plane, we 
attempt to solve the appropriate equation for a conserved probability density. 
It is also demanded that this quantity should agree with Born's probability 
rule along the real axis, which is the boundary condition in this case. An im- 
portant result obtained is that the characteristic curves of the conservation 
equation, along which the information about the solution propagates, are the 
same as the paths of particle in the complex trajectory representation. We find 
that there exist conserved probability densities which agree with the bound- 
ary condition in most of the examples considered. An alternative method to 
evaluate this in terms of a trajectory integral, in a manner similar to obtaining 
the probability distribution along the real line, is also presented. This 
latter method shall be of interest while solving the QHJE for the motion of 
individual particles. Most importantly, we show that the above conservation 
equation can be derived from this definition of probability, by using the com- 
plex extension of Schrodinger equation. 



2 Probability from velocity field 



Let us recall that in this scheme, the real part of the velocity of a particle on 
the real line, denoted 0), always agrees with the velocity of particle 

in the dBB representation pQ. To obtain the quantum probability function 
from the velocities, which is our first goal, we note further that the imaginary 
component of the velocity of the particle in the complex trajectory can be 
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written as 



h 

2m 



(5) 



This helps to write the probability density to find the particle around some 
point x = x r as 



2m 



ty*ty(x r , 0) = P(x r ) = A/"exp | — ^ J Xidx r J , (6) 



where the integral is taken along the real axis. This possibility of regaining 
the quantum probability distribution from the velocity field is a unique fea- 
ture of the complex trajectory formulation. For instance, in the dBB approach, 
the velocity fields for all bound eigenstates are zero everywhere and it is not 
possible to obtain a relation between velocity and probability. The FFM tra- 
jectory representation, on the other hand, does not claim any connection with 
probability. 

In the following, we limit ourselves to one dimension and to time-independent 
problems. The variable x is always assumed to be complex. As stated earlier, 
in this new representation, which is based on the existence of trajectories in 
the complex plane, even a particle in an eigenstate can be in any one of the 
infinitely many possible trajectories, depending on its initial position in the 
complex plane. Thus it is also desirable to look for the probability of a particle 
to be in an area dx r dxi around some point (x r ,Xi) in the complex plane. Let 
this quantity be denoted as p(x r , x.i)dx r dxi. 

An explicit expression for p(x r ,Xi) can be arrived at in the following way. 
Extending the probability density to the entire a; r a;j-plane demands a conser- 
vation equation of the form (for time-independent cases) 

d{px r ) + djpxj) = Q _ 
doc y> ^ 



To solve this partial differential equation, we may write p(x r , xi) = h(x r , Xi)f{p), 
where h(x r ,Xi) is some solution of Eq. (J7|) and p is some combination of x r 
and Xi, whose value remains a constant along its characteristic curves [19J. 
Substituting this form of p into (J7J), we see that the characteristic curves are 
obtained by integrating the equation 

dec if dx ij diXj'i cc a , . 

— = — or ;r^ = -, (8) 

Xj ip Xs^ (JjXj y Xj f 
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which is found to be the same as Eq. ([3]). This demonstrates the important 
property that the characteristic curves for the above conservation equation are 
identical to the complex paths of particles in the present quantum trajectory 
representation. 

We may now find the exact form of f(p) by requiring that p(x r , 0) agrees with 
the probability P(x r ), which is the boundary condition in this case. Let the 
integration constant in the above equation (jSJ) be the (real) coordinate of any 
one point of crossing of the trajectory on the real axis, denoted as x t q. Since 
the characteristic curves are identical to the complex paths, one can take x r o 
as the constant p along the characteristic curve and let it be expressed in terms 
of x r and Xj. The assumed form for the extended probability distribution p 
may then be written as 

p(x r ,Xi) = h(x r ,Xi)f(x rQ ). (9) 



Now we can choose f(x r o) subject to the boundary condition. At the point 
x = x r Q at which the curve C crosses the real line, we demand (the boundary 
condition) 

p(x r0 , 0) = h(x r0 , 0)f(x r0 ) = P(x r0 ) (10) 



and obtain /(x r0 ). Expressing x r o in terms of x r and Xi in /(x r0 ), Eq. flU]) gives 
p(x r , Xi). 

A word of caution is appropriate here. There may be instances, as we shall see 
below, when the boundary condition overdetermines the problem and we are 
unable to find a solution. It is observed that this happens in certain regions 
of the complex space where the trajectories were previously suspected to be 
nonviable. 

We may also note that the differential in Eq. ([3]), given by 

iidx r — x r dxi = (11) 



is inexact because 

dx^ dxj (12) 
doc y dx j 



unless both the partial derivatives are zero. This is due to the Cauchy-Riemann 
conditions satisfied by the analytic function x in the complex x-plane (except 
at its singular points). However, this inexact differential can always be made 
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exact by multiplying with an integrating factor p(x r , Xj), which obeys 
dx f doc <i 

Thus we see that the integrating factor p can serve as h(x r ,Xi) in Eq. (Ej). 

The following observation may be helpful in finding h(x r , X{), or even p(x r , Xi), 
in the case of some special potentials. Let us denote d^/dx = x( x )- Then using 
the time-independent Schrodinger equation, one can rewrite the equation of 
motion (J2J) as 

± = *£z!0x (14) 
h x 



where x ' = dx/dx. For V = Vq, a real constant, one can integrate this to 
obtain 



X{x) = A exp 



2i{E-V Q ) + 

h 



(15) 



from which we get the trajectory of the particle in the complex x-plane as 
X*( x )x( x ) H A | 2 - Moreover, one can see that 

(x0 y.iM$. (16) 

nr x*x 



But for constant potentials, (x ')*X ' °^ in the complex plane. It can be 
seen that the above expression f|T6|) satisfies the conservation equation and 
agrees with the boundary condition. In other words, the extended, conserved 
probability p(x r ,Xi) in this case can be written as 



p(x ry Xij \I/ \l/ . . , 

m z x*x 

and hence p varies inversely clS I X I 5 clS the particle moves along a particular 
trajectory with fixed A. This result is not contradictory to the WKB result that 
oc 1 /^classical f° r constant potentials [20], since | x \ along a trajectory 
in our case is not the same as ^classical' 

In the case of harmonic oscillator potential, we shall see below that an expres- 
sion of the form | A | 2 / (x*x) will give the solution h(x r , Xj) of the conservation 
equation, but for the particular solution which agrees with the boundary con- 
dition, one need to find f(x r0 ) too. 
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3 Examples 



Let us now first apply the procedure to the harmonic oscillator. In the n = 
case, we have i{ = (ha 2 /m)x r and the quantum probability density is regained 
by using Eq. (jBJ) as 



X 

2m 



P(x r ) = Afexp — ^ | x,rfx r = A/" e- Q2 < (17) 



To find p(x r , Xi), first let us note that h — ho, a constant, can be a solution of 
the conservation equation (j7|) in this case. Then we find f(x r o) using Eq. ffTUl) 
as proportional to exp(— a 2 £ 2 ). To put x r0 in terms of x r and Xj, Eq. ([3]) shall 
be integrated to obtain the paths as x 2 + x 2 = constant []]. Equating this 
constant to x 2 , we can finally write 

p(x r ,Xi) oce-" 5 ^. (18) 



This conserved probability is plotted in Fig. 1. 

Similarly in the n = 1 eigenstate of this case, one uses Eq. (J6]) to regain 
P(x r ) = M\x 2 exp(— a 2 x 2 ). In the next step, h(x r ,Xi) = [x 2 + xf) is found to 
be a solution of the conservation equation (One could guess this expression 
from the form of x*x in this it appears in its denominator.) f(x r o) 

can now be found as 

f(x r0 )^e- a2x ro. (19) 



The complex paths in this case are given by (a 2 x 2 — a 2 xf — l) 2 + 4a 4 x 2 x 2 = 
| A \ 2 , a constant. (Note that also this appears in the expression for x*x. 
Explicitly, one obtains x*x —\ A \ 2 j[x 2 + x 2 )). Equating this constant to 
(at 2 x 2 — l) 2 , one can obtain f(x r0 ) and also the extended probability density 
p, in terms of x r and Xi. But while taking square roots, one need to be careful. 
It shall be noted that for the region containing the subnests in the harmonic 
oscillator (with A < 1 in the present n — 1 case) PQ, the boundary condition 
overdetermines the problem, resulting in there being no solution. But for the 
region outside it, one can write 

f{x r0 ) oc exp(- v /(a 2 x 2 - a 2 x 2 - l) 2 + 4a*x 2 x 2 ), (20) 



and therefore, 

p(x r , Xi) oc (x 2 + x 2 ) exp(— J (a 2 x 2 — a 2 x 2 — l) 2 + 4a 4 x 2 xf ) (21) 
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which is plotted in Fig. 2. 



Now let us consider another example of a particle in an infinite square well 
potential, between x = and x = a. Using the energy eigenfuntion ip n = 
\j2/asm(mrx/a), one can find that the complex paths are given by 

coshf -J+cosf -J =| A | 2 . (22) 



We can choose h(x r ,Xi) = cosh(2mrxi/a) — cos(2mrx r /a) as a solution of the 
conservation equation. (Here again, we have x*x =\ A | 2 /h(x r ,Xi).) Interest- 
ingly, here, we have f(x r0 ) = constant and hence, 

p(x r ,Xi) (x cosh ( — cos ( . (23) 



The n = 1 case of this distribution is plotted in Fig. 3. 

A similar situation arises in the case of complex trajectories in the potential 
step problem discussed in pp. Again, we find x*x —\ A | 2 /h(x r ,Xi) and 
f(x r o) = constant. For the reflection constant equal to 1/2, we thus get, for 
the potential step, 

p(x r ,Xi) oc e~ 2kx * + >/2cos(2ifea; r ) + \e 2kx \ (24) 



which is plotted in Fig. 4. 



4 Trajectory integral method to find p 



It may now appear that the above two approaches to probability, namely, the 
one which gives the Born's \Er*\Ef probability along the real line using equa- 
tion and the other, which gives a conserved probability density over the 
entire avxj-plane by solving a conservation equation, do not have anything in 
common. But in this section, we show that there is a more elegant way of ob- 
taining p(x r ,Xi), by uniting them. This method may also be useful in finding 
the variation of probability along any trajectory. 

First, we postulate that if p , the extended probability density at some point 
(x r o,Xio) is given, then p(x r ,Xi) at another point that lies on the trajectory 
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which passes through (x r o,Xio), is 
p(x r ,Xi) = p exp 



^ J Im Qmx 2 + V(x) ) dt' 

to 



(25) 



Here, the integral is taken along the trajectory [x r (t'),Xi(t')}. One can de- 
rive the continuity equation by using the extended version of the Schrodinger 
equation, which gives 

( \ h dx\ 

Im(E) = Im -mx 2 + V(x) + = 0, (26) 

V 2 2iox 



since energy and time are assumed real 'J\ . This helps to write the above 
definition ( |25|) as 



dx r 



p{x r ,Xi) =p exp y 2 J ^~ dt 'j ' ( 27 ) 



which in turn gives 

dp _ _ 2 ^V^ _ _ ( + dxA 
dt dx r I dx r dx,i ) 



The last step follows from the analyticity of x. This leads to the continuity 
equation for the particle, as it moves along: i.e., 

dp + d(px r ) + d(pxj) = Q . 2g . 
dt dx r dxi 



While evaluating p with the help of (125j) above, one needs to know po at 
(x r0 ,x i0 ) and if we choose this point as (x r0 ,0), the point of crossing of the 
trajectory on the real line, then p may take the value P(x r o) and may be found 
using (jSJ). Here it shall be reminded that the integral in P(x r0 ) is evaluated 
over the real line only. 

To summarize the alternative method, we note that the conserved, extended 
probability density is 



p(x r , Xi) oc exp 



2m 

IT 



.t>0 



Xidx r exp 



— J Im Qmx 2 + V(x) ) dt' 
to 



(30) 
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with the integral in the first factor evaluated over the real line and that in 
the second factor over the trajectory of the particle. This method is easier to 
implement numerically on a computer. 

We can now demonstrate this trajectory integral method by applying it to the 
n = 1 harmonic oscillator state. In this case, taking the numerical values of 
both a and to to be unity (which then also applies to h/m), it is easily seen 
that pQ 

x r = ±\/(l + A cos 2t) ± vT + A 2 + 2A cos^/v 7 ^ (31) 



and 



J, ; 



A sin 2t 



(32) 



Here one can write P(x r o) = (1 + A) exp(— (1 + A)). Finding x r and Xi using 
the above expressions, we have evaluated numerically 

p(x r ,Xi) = P(x r0 )exp ( -4 / s x 2+ l x 2\2 dt I • 




for different A and t. This was found to be the same extended distribution as 
that given by Eq. (12TT) and shown in Fig. 2, with A > 1. More specifically, 
in this case, since equations ([9]) and (!30|) must agree, the surface plots of 
h(x r ,Xi) = {xi + x\) and the expression 



with various values of A, t but plotted against x r and Xi evaluated using 
equations (!3T!) and (|32|) . must be the same. This too is found to be true, 
which demonstrates that the two methods of evaluating the probability give 
identical results. 



5 Discussion 



The complex quantum trajectory representation is worth pursuing mainly for 
the alternative interpretation it offers to standard quantum mechanics. The 
standard quantum theory is spectacularly successful in explaining all obser- 
vations made so far and there is no demand for modifications to the quantum 
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probability axiom. It is envisaged that the new complex trajectory represen- 
tation does not differ from standard quantum mechanics, with regard to its 
predictions of experimental results. This is ensured by accepting the Born's 
probability axiom as such in the new theory, on the real line. The quantum 
probability density, though it is not claimed to be perfectly understood, works 
well in a miraculous manner. In this paper, we obtain this distribution along 
the real line from the velocity field of particles in the complex trajectory repre- 
sentation. Its conservation on the real line is guaranteed in standard quantum 
mechanics by virtue of the Schrodinger equation and there is no need for any 
separate proof for the same in the complex trajectory representation. What 
we try to do next is to extend it to the complex plane, where it should obey 
conservation laws, as in the case of any other probability function. Our great 
expectation is to find such a probability density in the complex plane, in the 
hope of it giving us a clue to the basis of this distribution. 

Here we have first adopted the strategy of solving the conservation equation 
in the x r Xj-plane, with the velocity field given by Eq. (jSJ), applicable to one 
dimensional, time-independent single particle states. In this endeavor, there 
is a significant improvement in our understanding of the origin of quantum 
probability. First, we have noticed that the characteristic curves of the conser- 
vation equation, along which the information about the solution propagates, 
are the same as the paths of particle in the complex trajectory representation. 
The most notable results we have obtained are that an extended, conserved 
probability p(x r ,Xi), which agrees with the quantum probability rule along 
the real line, can exist and that it can be derived from the velocity field of 
particles. There are parts of the plane in some examples where such a prob- 
ability cannot be found. But in those cases where it is possible, we find that 
p(x r , Xi) can be written as a product of two factors; one [f(x r o)], a constant for 
the given complex path and the other [h(x r ,Xi)], dependent on the velocity 
of the particle as it moves along this path. In particular, we note that the 
Born's probability rule = P(x r ) itself, which is valid on the real line, 
can be written as a product as stated above. On the other hand, in an alter- 
native method of evaluating the extended probability, p(x r ,Xi) is found to be 
obtainable ) times an exponential factor which involves an integral 

with respect to time along the trajectory, of the imaginary part of the sum of 
the complex kinetic and potential energies of the particle. The two methods of 
obtaining the extended probability density were compared and found to give 
identical results. The latter trajectory integral method is easier to implement 
using a computer and may be particularly useful while solving the QHJE for 
single particles. 

The solutions we obtain for the complex plane are not the same as except 
in the case of those constant potentials, discussed in Sec. (4). It shall be noted 
that p is defined only over that region of the x r Xi-plane, where characteristic 
curves which cross the real axis pass. At all other points, p is assumed to 
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be zero. The extended probability may be normalized on this basis. Also the 
present representation offers the viewpoint that a particle is observed only 
when it surfaces on the real axis. P(x r ) is the probability for such appearances 
and shall be used for computations involving measurable quantities. In this 
way, all the standard results in quantum mechanics continue to be unaffected. 

Most importantly, we have seen that the definition (125]) of the extended prob- 
ability density leads to the continuity equation (1291) . while using the complex 
extension of Schrodinger equation. This definition of p can thus be of general 
validity. But here, p is a real function of x r and Xj, which does not obey the 
Cauchy-Riemann equations, and hence is not analytic in the complex x-plane. 
It may be noted that in our attempt to have a better understanding of the 
origin of quantum probability rule, where we could make some inroads, this 
nonanalyticity does not appear as a problem. Ref. [JS] strongly advocates an 
analytic expression for p, in view of its anticipated advantages in the synthetic 
time-dependent Schrodinger equation applications. But the definition in [18] 
is shown not to lead to a continuity equation and it is presented there as a 
negative result. Whether our definition of extended probability in this paper, 
which does obey a continuity equation, is helpful for synthetic applications is 
an open problem, and is worth pursuing. 

To conclude, we summarize the positive and negative features of this possible 
definition of quantum probability density in the complex space, in compari- 
son with the distribution prescribed in [TSJ. First, we see that our definition 
helps to obtain the Born's probability density along the real line from 
the velocity field in the complex trajectory formalism, in a nontrivial way. 
On the other hand, the distribution in [TS] is simply the analytic continua- 
tion of and hence the agreement is trivial. Second, our definition of the 
probability density for the particle to be around some point in its trajectory, 
as it moves along, obeys a conservation equation. In contrast, it was shown 
in [T8J that the probability density defined in it is not conserved, in general. 
Furthermore, while solving the conservation equation, we have noted that p 
can be written as a product between two factors, one of which is a constant 
for a trajectory. This feature was deduced while attempting to solve the con- 
servation equation and is a consequence of the fact that trajectories and the 
characteristic curves of the conservation equation are identical in the complex 
trajectory representation. This has no parallels in the proposed probability 
density in [18]. A negative aspect of the present formalism is that the proba- 
bility density is not an analytic function in the complex plane, a feature which 
is suspected to adversely affect its applicability in the synthetic solutions. This 
is another feature that distinguishes the present distribution from the analytic 
one studied in [TS]. It is noted that in our attempt to understand the origin of 
quantum probability rule in complex trajectory representation, this nonana- 
lyticity does not pose any problem. But the complex probability that appears 
in [TS] off of the real axis would be an undesirable feature, and could have 
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posed some problems even if it was found to obey the continuity equation. 
Thus the issue of whether analyticity or real-valuedness of probability density 
is the more desirable feature in synthetic applications can be settled only when 
one progresses with such applications of the present formalism. 
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